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EQUATIONS  OP  THE  PROBLEM  OF  DETERMINING  THE  LOCATION  OF  A  MOVING 
OBJECT  3Y  GYROSCOPES  AND  ACC EL EEC METERS 


A.  Yu.  Ishlinsici*/  (Mcsccw) 

Summary 

The  problem  of  the  so-called  autonomous,  i.e.,  without  using 
external  references,  determination  of  the  location  of  a  moving  object 
is  of  great  practical  significance.  Until  recently,  this  problem  had 
virtually  no  chances  for  practical  resolution  because  of  the  lack  of 
precision  of  the  sensors  available  to  the  engineer  for  this  purpose, 
namely  -  newtonometers *,  or  accelerometers,  gyroscopes,  and 
integrators. 


Footnote:  |Re  are  proposing  that  the  instruments  previously  known  in 
technology  as  accelerometers,  or  acceleration  meters,  subsequently  be 
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called  newtonometers,  since  they  measure  the  combined  action  of  both 
the  forces  of  inertia  cf  translational  motion  and  the  forces  of 
gravity  on  their  sensor.  It  is  the  projection  of  the  resultant 
combination  of  these  forces  on  a  certain  direction  bound  to  the 
instrument,  wnich  we  will  subsequently  call  the  axis  of  sensitivity 
of  the  newtonometer,  which  is  measured. 

2e  will  point  out  tnat  the  above  forces  of  inertia  of 
translational  motion  shcula  be  determined,  of  course,  relative  to  a 
coordinate  system  which  is  strictly  bound  to  the  instrument  itself. 
Coriolis  forces  of  inertia  usually  do  not  affect  the  instrument 
reading.  End  footnote 


However,  recently  the  creation  of  new  elements  [1]  has  already  made 
it  possible  to  solve  this  prcolea  with  satisfactory  precision,  as 
long  as  the  duration  of  movement  of  the  object  is  short. 

This  report  considers  the  theory  of  one  possible  version  of  the 
problem  of  the  autonomous  determination  of  the  location  of  a  moving 
object. 


The  important  problem  of  estimating  the  precision  of  the 
determination  of  location,  whicn  is  affected  by  the  presence  of  the 
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so-called  instrument  errors  or  the  newtonometers,  gyroscopes  and 
other  conponents  of  the  system,  is  beyond  the  scope  of  this  report 
and  aay  become  the  object  of  special  studies.  Therefore,  in  the 
future  we  will  assume  that  the  above  components  operate  without 
errors.  Accordingly,  the  electromechanical  system  which  solves  the 
problem  will  be  considered  with  the  assumption  that  all  of  its 
parameters  precisely  correspond  to  their  theoretical  values,  and  that 
there  are  no  technological  errors  (a.g.,  lack  of  precision  of 
assembly  and  free  strokes  in  transmissions) .  Furthermore,  the  initial 
conditions  of  small  oscillations  or  this  system  can  be  arbitrary. 


1.  First  we  will  give  the  solution  cf  the  known  problem  of 
autonomous  determination  during  tne  movement  cf  an  object  over  an  arc 
cf  the  great  circle  of  a  certain  aonrotational  sphere  S,  whose  center 
coincides  with  the  center  of  tne  Earth  (Fig.  1)  .  In  the  simplest 
case,  this  corresponds  to  movement  at  a  constant  altitude  above  the 
equator.  Here  the  determination  of  the  location  of  the  object 
relative  to  the  Earth  itself  is  reduced  to  simply  timing  the 


movement 
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?ig.  1.  Fig.  2. 

We  will  bind  two  systems  of  coordinates  xy  and  t*7}*  having  the 
same  origin  with  a  certain  determined  point  of  a  sowing  object  (Fig. 

1)  .  Hanceforth,  we  will  caj.1  this  point  the  center  of  the  object.  Th8 
x-axis  of  coordinate  systea  xy  is  directed  along  the  velocity  vector 
of  the  object  v.  Accordingly,  tae  y-axis  is  the  continuation  of  the 
Earth's  radius  passing  through  the  center  of  the  moving  object.  The 
direction  of  the  axes  of  coordicate  systea  ?<“>£*  is  fixed  (relative  to 
stationary  stars  or,  analogously,  relative  to  sphere  S) . 

Suppose  that  the  soveaents  of  both  coordinate  systeas  coincide 
at  the  initial  point  in  time.  S uossquent ly ,  coordinate  system  xy  will 
be  turned  relative  to  progressively  moving  system  by  a  certain 


angle  *  which  varies  through  cine.  This  angle  is  related  to  the 
distance  s  =  s(t)  covered  by  the  center  of  the  object  froa  the 
initial  point  by  the  relationsmp 

?=7f  (1-D 

Here  H  is  the  radius  of  tae  arc  of  the  great  circle  (the  Earth's 
radius)  over  which  the  center  of  tae  object  moves. 

«Je  will  place  a  certain  platform  1  stabilized  by  gyroscopes  on 
the  object  (Fig.  2)  .  In  the  simplest  case,  this  platfora  is  oriented 
by  special  servo  systems  so  taat  it  is  always  perpendicular  to  the 
vector  of  the  intrinsic  kinetic  moment  of  a  certain  gyroscope  (Fig. 

2) .  The  bearings  of  the  axle  of  tae  outer  ring  of  cardan  joints  for 
the  suspension  of  this  gyroscope  can  be  attached  to  a  stabilizing 
platform. 

Let  the  intrinsic  Kinetic  moment  of  the  gyroscope  H  lie  in  the 
plane  and  a  certain  moment  h  =  !i(t)  be  applied  to  the  axis  of 

the  outer  ring  of  its  suspension.  In  this  case,  the  precession  of  the 
gyroscope  axis  occurs  in  the  same  plane  Following  the 

gyroscope,  the  stabilized  platform  will  rotate  at  an  angular  velocity 
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where  if/  is  the  angla  between  the  axis  of  natural  rotation  of  the 
gyroscope  and  axis 

Furthermore,  (Fig.  2)  we  will  place  nevtcnomster  2  in  the  plane 
of  the  stabilized  platform  1  and  we  will  uss  a  =  a(t)  to  designate 
its  current  reading. 

Let  the  axis  of  sensitivity  or  the  newtonometer  coincide  with 
the  x-axis.  Then  its  reading  ait)  (Fig.  3)  is  expressed  by  the 
formula 

a(t)  =  ^/—  —  Jain*  +j~£cos«  (13) 


Here  j  is  the  gravitational  force  per  unit  aass,  and  the  angle 
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«  =  tp  —  4»  (1.4) 

is  the  deflection  of  the  staoilizeu  platform  from  horizontal 
direction  (to  be  more  precise,  perpendicular  to  th9  Earth’s  radius). 
If  the  function  a(t)  obtained  at  tne  output  of  the  newtonometer  is 
integrated  using  special  devices,  as  a  result  we  can  fora  the  new 
function 

I 

K^a(t)dt+n,  (1.5) 

0 

where  K  and  m  are  certain  constant  parameters.  The  necessary  values 
of  these  parameters  will  oe  established  later. 

He  can  reproduce  function  (1.  5)  in  the  fora  of  moment  M  (t) 
acting  on  the  gyroscope  of  the  staoilized  platform.  Then,  plugging 
,1(t)  into  formula  (1.2)  and  integrating,  ws  will  have  ' 

II 

♦  =  7j  o(t)dt*  +  jjt  +  ^  (1.6) 

00 

where  g>0  is  the  angle  of  inclination  of  the  stabilized  platform  to 
the  horizon  at  the  initial  point  m  time  t  =  0  (to  be  more  precise, 
to  the  x-axis,  or  to  axis  €*,  whicn  coincides  with  it  at  this 
instant)  . 

Considering  formulae  (1.1)  anu  (1.4)  and  relationship  (1.6),  we 
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will  have  the  equation 

<i-7) 

00 

which  aakes  it  possible  to  sec  up  the  differential  equation  for  the 
function  a  =  a(t)  and  tne  initial  conditions  of  this  equation. 
Actually,  setting  t  =  0  in  equation  (1.7)  and  noting  that  (0)  =  0,  we 
will  have 

*  (0)  =  —  q*o  •  (i.8) 

which,  of  course,  also  rciiows  directly  from  formula  (1.4). 
Subsequently,  we  will  equate  tne  derivatives  of  the  left  and  right 
sides  of  equation  (1.7).  ee  will  octain 

tS-tSH “'9) 

0 

Whence  it  follows  that  tne  initial  value  of  the  derivative  of 
function  a  (t)  through  tine  is  expressed  by  the  formula 

dap  =  »Jp  _  *_  (, (0) -  (1-10) 

Here  v(0)  is  the  initial  value  of  the  velocity  of  the  center  of 
the  object  relative  to  ncnrotatmg  sphere  S. 


Finally,  if  we  differentiate  both  siles  of  equation  (1.7)  once 
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more  and  consider  formula  (1.3) ,  we  will  arrive  at  the  relationship 


3+v(/~x)— -(W«  .)£  d.»« 


With  a  given  function  s  =  s(t),  this  relationship  can  be 
considered  to  be  the  second-order  differential  equation  for  the 
function  a  =  at(t),  i.e.,  for  me  angle  of  inclination  of  the 
stabilized  platform  to  the  aocizos.  The  initial  conditions  of  this 
differential  equation  are  determined  by  equations  (1.3)  and  (1.10), 


2.  Differential  equation  vl.11)  has  a  partial  integral 


which  is  extremely  important  for  solving  the  problem  of  the 
autonomous  determination  of  the  position  of  a  moving  object,  as  long 


as  the  equation 


l  K_ 

“*  H 


is  satisfied,  as  it  is  net  hard  to  see,  and  the  initial  conditions 
are  such  that 


«(0)  -0, 


d2  (0)  n 

dt 


(2.3) 


c 
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0 


Equation  (2.2)  determines  the  value  of  parameter  K.  Together 
with  equation  (1.8),  the  first  condition  of  (2.3)  leads  to  the 
requirement  that  there  is  no  initial  deflection  of  the  stabilizing 
platform  from  the  horizon,  i.e.,  tne  equation 


^0  =  0 


is  satisfied. 


(2.4> 


According  to  the  second  condition  of  (2.3)  and  formula  (1.10), 
we  will  have 


^  (U1  _  HI 

H  ~  H 


(2.5) 


This  determines  the  value  or  parameter  a  in  the  device  which 
forms  function  (1.5). 

Thus,  when  conditions  (2.4)  and  (2.5)  are  observed,  as  well  as 
equation  (2.2),  the  stabilized  platform  will  remain  horizontal  (to  b<= 
■ore  precise,  it  will  he  perpendicular  to  the  Earth’s  radius)  with 
any  law  of  the  movement  of  the  center  of  the  object  s  =  s (t)  . 

However,  if  a  *  0,  a(t)  =  d*s/dtz  in  accordance  with  formula 
(1.3)  and,  consequently,  with  consideration  of  conditions  (2.5)  and 
(2.2),  function  (1.5)  will  assume  the  form: 
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A  +  »(°)]"^»(0  (2.6) 

0 

Thus,  in  this  case,  with  precision  down  to  the  constant  factor 
K,  function  (1,5)  is  the  currant  value  of  the  velocity  of  the  object 
v  *  v  (t)  relative  to  ncnrotatmg  sphere  S.  Therefore,  if  we  now  apply 
function  (1.5)  to  the  second  integrating  device,  with  the  same 
assumptions,  at  its  output  we  will  obtain  the  unknown  distance  s(t) 
covered  by  the  center  cr  tne  ooject  from  the  original  position  (also 
with  precision  down  to  tha  factor  ri)  . 

3.  In  tha  same  general  case,  when  a  fc.  0,  the  value 

<  i 

5  [&  \  a(t)dt  +  nAdt  (3.1) 

0  0  J 

generated  by  the  above  integrating  devices  differs  from  the  product 
Ks  (t)  .  According  to  equation  0*7)  and  condition  (2 •  2)  #  we  will  have 

i  t 

*  (0  —  iS  Ka  (l)  dl  ft*  (<)  +  ft*o  (3.2) 

0  0 

Thus,  the  value  of  tne  error  ns  in  determining  the  distance  by 


the  above  method  is 


As  -  «(<x(t)-*(°)l 


(3.3) 
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where  a(t)  is  determined  ny  dif xar antial  aquation  (1.11)  with  initial 
conditions  (1.3)  and  (1.10)  and  tne  observation  of  equation  (2.2). 

The  angle  <x  (t)  of  the  inclination  of  the  stabilized  platform  to 
the  x-axis  can  be  considered  to  be  extremely  small.  Therefore, 
dropping  the  terms  with  tne  second  order  of  saallnsss  relative  to  a 
from  equation  (1.11)  and  considering  equation  (2.2),  we  will  arrive 
at  the  following  homogeneous  linear  differential  equation: 

»+•*■('—£)—» 

With  the  known  approximation  we  can  assume  that 

_ _ —  s=  g  .=  const  (3-5) 

where  g  is  the  value  of  tae  acceleration  of  the  force  of  gravity  in 
the  region  of  movement  cf  tae  object.  In  this  case,  the  solution  to 
equation  (3.4)  is  in  the  fora  ox  a  geometric  function: 

«(f)  =  *(0)cosvt  +  ~  sin  v<  (v*  =>  -jjr)  (3.6) 

whose  period 

T  -  2k)/*77~«4.4  min. 
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is  known  in  gyroscope  tneory  *s  the  Schuler  period. 

Thus,  in  general,  according  to  formulae  (3.3)  and  (3.6),  the 
error  in  the  autonomous  determination  of  the  location  of  the  object 
by  the  above  method  fluctuates. 


4.  We  can  also  suggest  otner  aetaods  of  determining  the  location  of 
an  object  which  generally  lead  to  the  same  results. 

For  example,  suppose  the  platform  is  stabilized  relative  to 
coordinate  system  so  that  it  remains  parallel  to  the  €*-axis. 

This  stabilization  can  be  ootaiaed  by  free  gyroscopes  or  by  tracking 
stars  (astronavigation)  £2].  tais  case,  the  newtonometer  should  be 
turned  relative  to  the  platform  oy  the  following  angle 

*  =  +  (4.1) 

00 

for  which  its  readings  should  be  integrated  twice. 

With  the  precise  observation  of  certain  equations  similar  to  the 
preceding  equation,  we  nave 

«  =»  /?X  (4.2,^ 


There  are  original  devices  wnich  perform  the  doable  integration 
immediately  without  breaking  it  down  into  two  successive  operations 
(Boykov  integrator)  [  3  ]. 


5.  He  will  proceed  to  the  consideration  of  the  problem  of 
determining  the  location  of  an  ccject  whose  center  moves  randomly 
over  the  Earth's  sphere.  For  solving  this  problem  we  will  use  a 
platform  stabilized  by  gyroscopes,  so  that  the  angular  velocities 


around  the  x-  and  y-axes  lying  in  the  plane  of  the  platform  and  the 
z-axis  normal  to  it  arose  oecausa  three  moments  Ht,  n2  and  applied 
to  the  axes  of  the  gyroscope  housings. 


Pigure  4  shows  a  possiole  diagram  for  realizing  this  type  of 
stabilization.  The  plane  of  the  platform  ?,  which  is  connected  with 
the  object  by  a  three-dimensional  cardan  suspension  (not  shown  in  the 
figure),  is  continuously  made  to  coincide  with  the  plane 
perpendicular  to  the  intrinsic  axis  of  the  gyroscope  1  by  servo 
systems.  The  bearings  cf  the  axis  of  the  outer  ring  of  this 
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gyroscope's  cardan  suspension  are  arranged  on  the  stabilized  platfora 
so  that  the  axis  itself  lies  in  tne  plane  of  the  platfora.  When  the 
servo  system  operates  perfactiy,  the  axis  of  the  inner  ring,  i.e., 
the  axis  of  the  gyroscope  nousing,  also  coincides  with  the  plane  of 
the  platfora.  Let  fSt  and  a2  be  tne  moments  applied  to  the  axis  of  the 
outer  ring  of  the  suspension  and  tne  axis  of  the  housing, 
respectively,  designated  ay  x  and  y.  The  moments  and  Ji2  cause  the 
precession  of  the  gyroscope  and,  consequently,  the  rotation  of  the 
platfora  around  the  x-  and  y-axes  at  angular  velocities  and 
The  values  of  the  latter  ars  deterained  by  the  first  two  formulae  in 
(5.  1)  . 


Pig.  4. 


The  bearings  of  the  axis  of  tne  outer  cardan  ring  of  the 
gyroscope  2  (Pig.  4)  are  connected  to  the  platfora  with  the  saee 
intrinsic  kinetic  aoaent  d  as  tne  first.  The  axis  of  this  outer  ring 
is  perpendicular  to  the  plane  of  tne  platfora.  A  certain  correcting 
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moment  is  applied  to  this  axis  so  that  the  natural  rotation  axis  of 
the  gyroscope  1  is  parallel  to  the  plane  of  the  platform. 
Furthermore,  a  special  servo  system  continuously  makes  the  above 
x-axis,  which  should  be  considered  to  be  rigidly  bound  to  the 
stabilized  platform,  coincide  with  the  natural  rotation  axis  of  the 
gyroscope  2.  The  moment  applied  to  the  axis  of  the  housing  of  the 
second  gyroscope  causes  tie  precession  of  the  gyroscope  and, 
consequently,  the  rotation  of  tae  platform  parallel  to  the  z-axis.  at 
angular  velocity  uz.  The  latter  is  perpendicular  to  the  plane  of  the 
platform,  and  together  with  tne  x-  and  y-axes,  forms  rectangular 
coordinate  system  xyz,  whicn  is  rigidly  bound  to  the  platform,  tfe 
will  consider  the  origin  or  tnis  system  to  be  located  in  the  center 
of  the  object.  The  values  of  ti3  and  are  related  by  the  third 
formula  of  (5.1)  . 

Two  newtonometers  whose  axes  of  sensitivity  are  directed  along 
the  axes  of  coordinates  x  and  y  are  located  on  the  stabilized 
platform  in  the  direct  proximity  or  the  origin  of  coordinate  system 
xyz.  Let  the  center  of  the  moving  object  move  randomly  over  the 
Earth's  sphere,  and  the  moments  and  n2  related  to  the  readings 
of  the  corresponding  newtonometers  and  a^  by  the  formulae 

I  ( 

Mt  =  —  K  —  m,t  Mt  u-s  K  ^  a^dt  -f  m,  (5.2) 

0  9 

for  which  special  integracing  devices  must  be  provided. 


I 
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We  will  explain  what  parameters  K ,  at  and  m2 ,  and  also  moment  n3 
aust  be  in  order  for  the  plane  or  the  platfora  to  remain  horizontal. 


Tha  projections  of  tae  acceleration  of  the  origin  of  coordinates 
relative  to  nonrotating  sphera  S  on  the  x-,  y-  and  z-axes  is 
expressed  by  the  formulae  [  4,  5] 

dv 

Wx  —  (-  UySl  OifVy 

iVy  =  +  oM>*  —  oixvt  (5-3) 

dvz 

W  j  =  -  -  4“  WjfBy  Uyl>x 

where  vx,  and  v^  are  t^e  projections  of  the  velocity  of  the  origin 
of  coordinates  on  the  sane  axes  relative  to  sphere  S.  In  the  case  in 
question,  vx  *  0,  and  the  rorca  cf  gravity  has  projections  equal  to 
zero  on  the  x-  and  y-axes.  Therefore,  the  newtono«eters  directly 
measure  accelerations  Wj,  and  w^,  and  according  to  foraulae  (5,1)  and 
(5.2),  the  following  equations  hold: 


o 

**  =  ““•»«)  *  +  Tr 


(5.4) 


Since  the  platfora  aust  remain  horizontal  and  the  z-axis 
directed  along  the  Eartn*s  raaius  as  the  object  aoves,  in  aquations 
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(5.4)  we  should  [5]  consider  mat 

VX  =  OtyR,  Vy  —  —  (J>XR  (5.5) 


Using  these  formulae,  we  wall  eliminate  the  values  wx  and 
from  equations  (5.4) .  He  will  nave 


vx 


KR  r 

(  dvx 

H  \ 

V  dt 

• 

**  r( 

(d"v 

t ry 

\  dt 

0 


—  vv)dt  -j-  -ji  m* 
+  dt  -j-  ~  m1 


(5.6) 


aquations  (5.6)  must  ne  satisriad  during  a  random  change  in  vx 
and  v^,,  i.  e. ,  they  must  be  identities.  However,  this  is  only  possible 
with  the  following  conditions: 

^L=l,  w,  =  0,  mj  =  A'rw(0).  m,  =  Kvx  (0)  (5.7) 

n 


The  first  condition  coincides  with  equation  (2.2);  the  second, 
according  to  the  third  formula  of  (5.1),  leads  to  the  requirement 

V,  =  0  (5.8) 

Finally,  the  last  two  conditions  concern  the  agreement  of  the 
initial  velocity  of  the  center  or  the  moving  object  (relative  to 
sphere  S)  ,  or  to  be  more  precise,  its  projection  on  the  x-  and 
y-axas,  with  parameters  nt  and  m2  of  the  integrating  devices. 
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6.  According  to  formulae  (5.1),  (5.2)  and  (5.3),  as  a  result  of 

the  work  of  the  integrators,  tad  projections  of  the  angular  velocity 
of  the  stabilized  platform  »»x  and  are  now  known  time  functions, 

=  0,  and  the  z-axis  is  directed  along  the  Earth's  radius. 
Theoretically,  this  makes  it  possiole  to  continuously  determine  the 
location  of  the  object  on  zaa  harm's  sphere  and  its  course. 

Actually,  we  will  consider  (Fig.  5)  the  so-called  geographic 
triangle  (t,;,  whose  apex  coincides  with  the  origin  of  coordinate 
system  xyz  (i.e.,  with  the  center  of  the  object);  side  5  is  directed 
towards  the  east,  side  ^  -  towards  the  north  and,  finally,  side  £  - 
from  the  Earth's  radius  upward.  The  projections  of  the  angular 
velocity  u  of  this  triangle  relative  to  nonrotating  sphere  S  on  sides 
€,  7}  and  £  are  expressed  £6]  by  taa  formulae 

y  y 

“i*3  ^  =  —jf  +  U  cos  <p,  u-  **  -jdL  tg  tp  -f-  U  sin  <p  (6.J£) 
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Fig.  5. 


N 


Fig.  6. 


Here  0  is  the  angular  velocity  of  the  Earth,  *  is  the  current 
latitude  of  the  location,  ana  1B  and  are  the  eastern  and  northern 
components  of  the  velocity  of  tne  center  of  the  object  relative  to 
the  Earth,  respectively.  Ooviously, 

VE-Jtcoa9£,  VN  =  Jl£  (6.2) 

where  X  is  the  current  longitude. 

In  this  case,  side  £  coincides  with  the  z-axis.  Therefore,  the 
projections  wy,  and  are  related  to  the  values  and  by 

the  relationships 


L  * 
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w*  =*  cos  x  4-  sin  x 

<i»v  =  —  Uj  sin  x  +  cos  X>  “» *■  “t  +  (®-®) 

where  X  is  the  angle  between  tna  €-  and  x-axes  (Pig.  6). 

Considering  formulae  (o.  1)  ana  (6.2)  and  the  fact  that  =  0, 
and  that  «*  *»  <^*(0 .  and  uw  =««(*)  are  unknown  time  functions,  according  to 
relationships  (6.3)  we  will  arrive  at  a  system  of  three  differential 
equations: 

—  cosx  +  (U  +  -^)cos<p  sinx==«x(0 
$8inX  +  (U  +  cos?  cos  x»«„(<)  (6.4) 

(£/+$)sin? 

for  the  three  unknown  functions  w(t),  X(t)  and  X  (t)  . 

With  known  initial  conditions  0(0)  ,  X(0)  and  X(0),  i.e.,  data  on 
the  location  of  the  object  and  its  orientation  at  the  initial  point 
in  time,  equations  (6.4)  can  oe  integrated  on  a  special  computer. 

Here  it  is  helpful  to  solve  taea  first  for  the  derivatives,  i.e.,  to 
represent  then  in  the  fora 
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=  —  o),  (/)  cos  x  +  ««(0  sin  X 

A  _  ,,  ,  «■»,(«)  »iPX  +  “y(0«»X  (6.5) 

'dt  ~~  ^  cos  9 

=  -  («,  (i)  sin  X  +  “y  (0  cos  X)  ® 

af 

Having  determined  functions  *  (t)  and  x  (t)  ,  we  can  also  find  th 
course  of  the  object,  i.e.,  the  angle*  which  its  velocity  vector 
foras  relative  to  the  Earth's  surface  with  the  meridian  of  the 
location  (Pig.  7).  According  to  formulas  (6.2)  and  (6.5),  the  cours 
fc  is  determined  by  the  equation 

U  co s  <P  —  {<)  sin  X  —  (1)  COS  x 

t  (7  ti  =s  - — — : -  (0.0 

K  wx  (1)  cos  x  —  «•»„  (i)  Sin  X 


iy  ar£h 


Fig,  7 . 


7.  now  we  will  proceed  to  studying  the  saall  movements  of  the 


I 
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stabilized  platform  with  tne  assumption  that  at  the  initial  instant, 
its  plane  was  not  in  a  precisely  horizontal  position  and  conditions 
(5.7)  of  the  selection,  of  tae  value  of  parameters  mt  and  m2  have  been 
satisfied  with  a  small  error.  Purmermore,  we  will  consider  that  the 
first  two  equations  of  (5.7)  nave  ceen  realized  with  complete 
precision. 


We  will  introduce  the  Darooux  triangle  x°y°z0,  whose  sides  x° 
and  y°  are  tangent  to  tae  Kama's  sphere  [5,  7]  and,  consequently, 
also  nonrotating  sphere  S.  Side  x°  is  directed  along  the  velocity 
vector  v  of  the  apex  of  tae  triangle  relative  to  sphera  S.  We  will 
call  triangle  x°y°z°  a  natural  Darnoux  triangle,  we  will  place  its 
apex  in  the  center  of  the  moving  object,  i.e.,  at  the  origin  of 
coordinate  system  xyz,  whica  is  rigidly  bound  to  the  stabilized 
platform.  The  projections  of  tae  angular  velocity  o°  of  the  natural 
Darboux  triangle  on  sides  x®,  y°  and  z°  are  represented  by  the 
formulae 

i^.0  =  0,  Wy.°  = -f- ,  w,.°  =  w  (7.1) 


At  a  given  velocity  v  =  v(t),  functions  w  =  S(t)  determine  [7] 
the  geodesic  curvature  cf  the  trajectory  of  the  apex  of  triangle 
x°y°z®  on  sphere  S. 


L 


DOC  =  0925 


PAG£  24 


In  turn,  the  projections  or  tne  acceleration  of  the  apex  of  the 
triangle  during  its  movement  relative  to  nonrotating  sphere  S  onto 
sides  x°,  y°  and  z°  are  expressed  by  the  formulae 

U)z-  =  ~~ ,  wv-  =Zv,  wz*  —  (7.2) 

Henceforth,  we  will  use  tne  table  of  cosines  of  the  angles 
between  the  axes  of  coordinate  system  xyz  and  the  sides  of  the 
triangle  x°y°z0.  It  is: 

x° 

a  *0*  T  CO(  3  —  SiU  Y  tlD  a  SlD  @ 

•  —  Y  CO*  a 

I  OO*  Y  »IO  @  +  *10  Y  »to  «  CO*  fl 


•to  Y  CO*  B  +  CO*  Y  ItD  1  *10  3 
CO*  Y  CO*  a 

*!n  y  Ho  B  —  oo*  y  »lo  »  co*  B 


—  co*  a  tin  B 
•to  * 

co*  a  oo*  B 

(7.3) 


Here  the  angles  a,  7  describe  the  orientation  of  coordinate 
system  xyz  relative  to  triangle  x°y°z°  (Fig.  8) .  The  angle  y  is  the 
angle  of  rotation  of  the  auxiliary  coordinate  system  x’y’z*  relative 
to  the  triangle.  The  z'-axis  of  this  system  coincides  with  side  z°. 
The  coordinate  system  is  rotated  counterclockwise  (viewed  from  the 
positive  part  of  the  z'-axis)  until  the  x*-axis  coincides  with  plane 
zx.  Similarly,  angle  a  is  determined  by  the  relative  position  of 
coordinate  system  x’y’z*  ana  another  auxiliary  system  xny,*zn,  the  x'- 
and  x"-axes  of  whica  coincide.  Tne  z"-axis  of  the  latter  coordinate 
system  also  coincides  wita  plane  zx,  as  a  result  of  which  its  axis  y" 
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turns  out  to  be  directed  along  toe  y-axis. 


Fig.  d. 


Whan  a  >  0,  coordinate  system  x*y"z*  is  turned  counterclockwise 
relative  to  system  x'y'z',  if  we  observe  its  rotation  from  the 
positive  part  of  the  x '-axis  (or,  analogously,  the  x“-axis) . 

Finally,  angle  0  is  tna  angle  between  x-  and  x*-axes  of 
coordinate  systems  xyz  and  z»y»z«t  The  sign  of  angle  ^  is  determined 
analogously  to  the  signs  of  angles  y  and  a. 

The  angular  velocity  w  of  coordinate  system  xyz  relative  to 
nonrotating  sphere  s  is  the  geometric  sum  of  the  angular  velocity  w° 
of  the  natural  Darboux  triangle  relative  to  the  same  sphere  and  the 
three  relative  angular  velocities  dy/dt,  da/dt  and  d0/dt.  The  latter 
are  the  angular  velocities  of  coordinate  system  x'y'z'  relative  to 


.A. 
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triangle  x°y°z0,  of  system  xrty,,zl*  relative  to  x'y'z',  and,  finally, 
of  coordinate  system  xyz,  i.e.,  tha  stabilized  platform,  relative  to 
system  x"y"z". 

The  relative  angular  velocity  dy/dt  is  diracted  along  side  z°, 
the  angular  velocity  da/dt  -  alcng  the  x'-axis,  and  angular  velocity 
d3/dt  -  along  the  y-axis.  The  x'-axis  coincides  with  side  x°  with  an 
angle  y  =  0  which,  according  to  Taale  (7.3),  makes  it  possible  to 
find  the  cosines  of  the  angles  formed  by  the  direction  of  the 
relative  angular  velocity  da/dt  to  the  axis  of  coordinate  system  xyz 
(Fig.  8)  .  Considering  all  of  the  acove,  we  obtain  the  following 
expressions  for  the  prelections  ux,  u^and  wfcof  tae  angular  velocity 
of  the  stabilized  platform  onto  the  axes  of  coordinate  system  xyz, 
which  is  rigidly  bound  tc  it: 

<0Z  =  -g-  (sin  y  cos  p  -f  cosy  sin  a  sin  p)  -f  (—  cos  a  sin  P)  -f-  ^  cosp 

cosy  cos  a  +  (w  +  <£')  sin  a  +  ^  (7.4) 

<o,  =  -jp  (sin  y  sin  p  —  cosy  sin  a  cos  P)  (u>  +  j  cos  a  cos  p  sinp 

Using  table  (7.3)  and  formulae  (7.2),  we  then  arrive  at  the 
expressions  for  the  projections  wx  ,  w^  and  of  the  acceleration  of 
the  origin  of  coordinate  system  xyz,  namely 
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(cosycosp  —  siiiysinxsin  p)  -f-  t*w>(sinycosp  +  cosysinxsin P)  — 

jjk 

—  cos  x  sin  p) 

(—  sinycosx)  -f- ow  cos  y  cos  x  —  ^-sinx  (7.5) 

(cos  y  sin  p  +  sin  y  sin  a  cos  P)  -j-  «3»  (sin  y  sin  p  —  cos  y  sin  *  cos  P)  — 

V*  Q 

—  cos  x  cos  p 

In  this  casa,  unlixe  50,  tna  projections  of  the  force  of  gravity 
on  the  x-  and  y-axss  are  aireaay  nonzero.  Therefore,  the 
nawtonometers  located  cn  tne  x-  ana  y-axes  will  measure  the  sum  of 
accelerations  w^  and  corresponding  to  the  projections  (with  the 
opposite  signs)  of  the  acceleration  of  the  force  of  gravity  j  onto 
these  same  axes.  This  acceleration  is  directed  along  side  z°  toward 
the  canter  of  the  Earth.  Now,  considering  formulae  (7.5)  and  (7.3), 
we  will  find  that  the  readings  ax  and  a^  of  the  corresponding 
newtonometers  should  be  expressed  ny  the  formulae 

<i*  =  ^  (cos  y  cos  p  —  sin  y  sin  xsin  p)  +  <2i>  (sin  y  cos  p  -f  cosy  sin  a  sin  P)  + 

+  (/— gy)(— cos«sinP)  (7.6) 

av  =  —  (—  sin  y  cos  a)  +  cos  y  cos  x  +  (j  —  jfj  sin  x 

Now  we  will  form  tne  moments  ilt  and  n2  according  to  formulae 
(5.2)  ,  while  n3,  like  before,  will  be  considered  to  be  equal  to  zero. 
According  to  relationships  (5.1),  we  will  have 


dv 

W  x  =  -r~ 

1  dt 


wu  = 


W,  - 


dv 

dt 

dv 

nr 


j 
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,  <  i 

“« ■  ~  4  5  -  -Jr  •  “«  -  -jr  j  axJi  +  ,  «,  =  o  (7.7) 

•  o 

As  follows  from  formulae  (7.4)  and  (7.6),  these  equations  can  be 
considered  to  be  the  system  of  equations  for  determining  the  time 
functions  a  (t)  ,  j3(t)  and  yit)  wita  the  given  functions  v(t)  and  «(t). 

Differentiating  the  first  two  equations  of  (7.7)  with  respect  to 
time,  we  will  arrive  at  a  system  of  differential  equations  for  the 
same  unknown  functions: 

^  +  -5-.-°.  w‘~°  <7-8> 

Considering  angles  a  and  0  in  formulae  (7.4)  and  (7.6)  to  be 
small,  and  disregarding  tae  second-order  terms  relative  to  these 
variables,  we  will  reduce  cae  expressions  for  u*, <*>»,<*>», <»x  and  to 

«*- K-wnY-^+aijP  +  ar 

v  ,  ( .  dy  \  .  rffS 

«w  =  s-COST  +  (m  +  -I)ci  +  J 

<*>z  —  ^(Psiny  —  «cosy) +«-+■  ^  (7.9) 

a*  =  jfcosr  +  ,:‘v  8iur—  (/“3f)p 

v°v  “  -  57  sin y  4-  <'»  cosy  +  (/  -  £)« 


Since  2.  0,  according  to  tne  third  formula  of  (7.9), 


DOC  *  0925 


PAGE  29 


‘-  +  17  =  S-(«cosY-M‘i»Y)  (7-10) 

Whance  it  follows  that  taa  terms 

(“+$)«.  (-+£>  '  <7-»> 

contained  in  the  first  two  formulae  of  (7.9)  are  of  the  second  order 
of  smallness  relative  tc  verianles  <r  and  3  and,  consequently,  they 
can  be  Iropped.  Thus,  wita  precision  down  to  the  second  order  of 
smallness,  we  will  have 

V  da  V  .  d  ^  (7  inv 

=  S'“Y  +  17.  COSY +  ^7  (7-12) 

Mow,  if  w e  use  relationship  (7.10)  to  eliminate  function  3(t) 
from  the  fourth  and  fifth  formulae  of  (7.9),  after  simple 
transformations  we  will  ootain  the  following  expressions  for  the 
values  of  a^  and  a^: 

«»“ar(*® 0SY)+  ^  (« cosT sin Y) sin Y“  (/-■?■)?  ?  ^  . 

av  =  —  ~  (»  sin  y)  +  -jy  («  cos  y  —  P  sin  y)  cos  Y  +  (/  y  )  * 

Finally,  substituting  tne  values  of  ax  and  a^  in  the  first 

two  equations  of  (7.8)  according  tc  formulae  (7.12)  and  (7.13)  and 
considering  equation  (2.2),  after  obvious  simplifications,  we  will 
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arrive  at  two  second-order  differential  equations: 

~  +  Jf  *  =  £  (*sio^  +  P  C0SY)  sin  Y 
g  +  =  ^l(«sinY+  ?  cosy)  cosy 


Together  with  relationship  (7.10),  they  for*  a  systea  of 
dif farantial  equations  fcr  tna  functions  a(t),  0(t)  and  y[ t)  . 

Considering  tha  saalinass  of  the  values  cf  a  and  0,  angle  y  in 
equations  (7.  14)  can  be  replaced  Dy  the  integral 

I 

Y  *  —  ^C'dt  -f-  Yo  (7.15) 

t 

where  y0  is  the  value  cf  tne  angle  y  between  the  x-axis  and  side  x° 
(i.e. ,  the  velocity  vector  of  tae  center  of  the  object)  at  the 
initial  point  in  tine.  As  it  is  easy  to  confirm  from  equations  (7.  10) 
and  (7.14),  the  errors  which  arise  during  this  substitution  are  of 
the  second  order  of  smallness. 

Thus,  the  study  of  small  oscillations  of  a  stabilized  platfora 
is  reduced  to  the  integration  of  two  linear  differential  equations 
with  variable  coefficients  depending  on  time1. 
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Pootnote:  *It  is  curious  to  note  tae  complete  identity  of  equations 
(7.14)  of  small  oscillations  of  the  stabilized  platform  in  question 
with  the  equations  of  small  oscillations  of  a  certain  physical 
pendulum  vith  a  reference  point  moving  over  sphere  S.  The  equilibrium 
conditions  of  this  pendulum  relative  to  the  natural  Darboux  triangle 
are  given  in  report  [5],  End  footnote 


8.  Setting  aside  the  integration  of  the  differential  equations 
(7.14)  of  small  oscillations  of  tae  platform  during  random  movement 
of  the  center  of  the  object  ever  tne  Earth's  sphere,  we  will  limit 
ourselves  to  the  case  u  -  u  and  v  =  const,  which  corresponds  to 
movement  at  a  constant  velocity  over  an  arc  of  the  great  circle  of 
nonrotating  sphere  S.  In  tais  case,  movement  relative  to  the  Earth 
will  occur  with  a  variable  relative  velocity  over  a  complex 
trajectory.  Setting  «  =»  0  in  formula  (7.  15)  ,  we  conclude  that  the 
angle  y  in  system  of  differential  equations  (7.14)  should  be 
considered  to  be  constant.  As  a  result,  system  (7.14)  is  broken  down 
into  two  independent  equations: 

£  (.  »» •,  + 1>  «>• r>  +  (i  -  i )  (* 1 t  ' +  t  T)  -  0  (8  () 

(«  cos  y  —  P  sin  y)  +  -jf  (*  cos  Y  P  s*n  Y)  "*  ® 

The  first  of  them  corresponds  to  the  angular  oscillations  of  the 
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platform  around  the  side  y°  of  taa  natural  Darboux  triangle,  and  the 
second  -  around  side  x°,  respectively .  As  we  already  mentioned  above, 
side  x°  is  directed  along  velocity  vector  v.  The  frequencies  of  these 
oscillations  are  close  together,  as  long  as  the  value  of  the  velocity 
v  is  not  too  great  (e. g,  ,  a  value  less  than  the  velocity  of  points  on 
the  Earth's  equator  during  its  diurnal  rotation).  This  corresponds  to 
a  period  of  time  approxixately  equal  to  84  min.  (i.e.,  the  Schuler 
period) . 


1 


9.  In  the  presence  of  small  oscillations  of  a  stabilized 
platform,  i.e.,  vhen  angles  a  and  p  are  nonzero,  relationships  (5.5) 
will  be  realized  with  a  certain  error.  Furthermore,  the  initial 
equations  (5.6)  for  this  method  of  autonomous  determination  of 
location  cannot  be  turned  into  precise  identities  because  of  the 
presence  of  projections  of  tha  acceleration  of  gravity  j  in  the 
newtonometer  readings  a^  and  a^.  one  would  expect  the  errors  in 
determining  the  latitude  and  longitude  of  the  location  of  the  moving 
object  and  its  course,  which  occur  because  of  the  above 
circumstances,  to  fluctuate.  However,  additional  research  is 
necessary  for  the  precise  determination  of  the  nature  of  the  change 
in  these  errors  through  time. 
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10.  Above  it  was  assumed  that  the  center  of  the  object  lores 
over  the  sphere  S,  on  tae  basis  of  which  it  was  assumed  that  v^  -  0 
inforaulae  (5.3).  Now  we  will  shew  how  to  eliainata  this  restrictive 
condition. 


Suppose  that  a  platxors  staoilized  by  gyroscopes  aoves  so  that 
its  plana  reaains  perpendicular  to  the  Earth's  radius.  In  this  case, 
it  is  necessary  to  satisfy  the  same  equations  (5.5),  where  v*  and  v^., 
li)ce  before,  are  the  projections  ox  the  velocity  v  of  the  center  of 
the  cardan  suspension  or  zaa  piatfera  onto  the  x-  and  y-axes  relative 
to  the  sphere  S,  and  t*A  and  ar*  the  projections  of  its  angular 
velocity  u  onto  the  saae  axes.  Thus,  we  should  have 


«x  * 


(10.1) 


(Jnlike  equations  (5.5),  aere  S  *  a  (t)  is  the  variable  value  of 
the  distance  between  the  center  of  the  cardan  suspension  of  the 
platform  and  the  center  cf  spaere  S. 

The  angular  velocities  wK  and  w^.  are  created  by  applying  aoaents 
nt  and  to  the  platfcra  according  to  the  first  two  foraulae  of 
(5.1).  Consequently,  ocaeats  and  H2  should  be  focaed  according  to 
the  aquations 
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Mx  —  — VV‘  =Tiv‘  (10  2) 

We  will  place  two  newtonometers,  whose  axes  of  sensitivity  are 
directed  along  the  x-  and  y-axas,  which  lie  in  the  plane  of  the 
platfora,  in  direct  proximity  to  tne  center  of  the  cardan  suspension. 
Their  readings  a^ft)  and  a^,(t)  will  not  contain  the  projections  of 
the  acceleration  of  gravity  j  onto  the  x-  and  y-axas,  since  by 
assumption,  as  it  moves  taa  stabilized  platform  remains  perpendicular 
to  the  straight  line  connecting  the  center  of  the  cardan  suspension 
and  the  canter  of  sphere  3.  Based  on  this,  we  can  set 

ax  (t)  =:*>*=  —  +  <**»*-  «*»..  “V  ('1=  “'v  =  dT  +  -  a*v  «  (10-3> 

Hare,  unlike  the  case  in  96,  vx  is  already  nonzero  and  is 
expressed  as 

_  (10.4) 

p*  •  dl 


Like  before,  if  we  consider  Hj  =  0,  then  according  to  the  third 
formula  of  (5.1),  we  again  obtain  0,  i.e«,  the  platfora  will  not 

have  its  angular  velocity  component  u  along  the  z-axis  directed  along 
the  radius  of  sphere  S.  Considering  this,  as  well  as  foraulae  (10.1) 
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and  (10.4),  according  tc  aquations  (10.3),  va  arrive  at  the 
relationships 


flx(J)-  Jr  +  T~, 


,  m  =  ,  •’v  dR 

*  +  7T3T 


(10.5) 


These  relationships  can  ae  considered  to  be  differential 
equations  which  make  it  possiole  to  plot  functions  vA  (t)  and  v^.(t) 
with  known  functions  a^,  (t)  and  a^(t)  ,  as  well  as  a  *  R  (t)  .  The  values 
of  the  former  functions  are  necessary  for  forming  moments  and  M2, 
which  control  the  orientation  of  the  platform. 

The  solution  of  equations  (10.5)  is  reduced  to  the  quadratic 
eq  nations 

i  i 

V'  =  ±  [\  Ha  (i)  dl+B  (0)  v,  (0)].  vv  =»  -i-  [[Bay  ( t )  dt+R{ 0)vy  (0)](10.6) 

«  o 


Thus,  according  tc  formulae  (10.2),  moments  (1,  and  M2  should  be 
formed  from  the  nevtoncmeter  readings  aA(t)  and  a^(t)  as  follows: 


Mt  * 
M v 


"  it  G R  w  w  dt + R  (0)  ry  (0)1 

^0 

[J  Jl  (t)  a,  (0  dt  +  R  (0)  r,  (0)] 


(10.7). 


For  this  purpose,  devices  which  multiply  and  divide  the  current 
values  must  be  provided,  as  well  as  integrators. 
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The  variable  value  cr  d  =  a  (t)  in  equations  (10.7)  is  considered 
be  known  in  advance.  However,  tne  presence  of  a  third  newtonometer 
whose  axis  of  sensitivity  is  parallel  to  the  z-axis  (the  Earth's 
radius)  theoretically  makes  it  possible  to  determine  this  function 
independent  of  any  other  auxiliary  devices.  Actually,  the  reading  of 
this  newtonometer  is  determined  by  the  formula 


a,  «  vx  -  /=  -5p  +  <-.4  -  0ivVx  -  1 


(10.8) 


Here 


,  V 

l)o 


(10.9) 


the  acceleration  of  gravity,  wnica  decreases  with  the  increase  in  the 
distance  from  the  Earth's  canter,  j0  is  its  value  on  the  Earth's 
surface,  whose  radius  is  designated  as  30. 


Using  formulae  (10.1),  (10.4)  and  (10.9),  we  obtain  the 

following  for  function  a  (t) 

Z*r  •*'  +  *,*  a  (t)  (10.10) 

- - -  Ht  atKi) 


The  device  which  integrates  tnis  differential  equation  must  be 
connected  to  on?  system  witn  integrators  which  reproduce  moments  fi2 
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and  .1 , .  According  to  formulae  (10.2),  those  moments  differ  from 
functions  v^(t)  and  v,  (t)  only  in  the  constant  factor.  Actually, 
functions  Y^(t)  and  v^(t)  are  m  equation  (10.10);  in  turn,  function 
R(t)  is  used  in  formulae  (10.  o). 


Problems  of  the  staciiity  cr  this  computational  system  require 
special  examination,  without  discussing  the  study  of  small 
oscillations  of  the  stabilizing  platform,  either,  we  will  point  out 
that  the  subsequent  solution  or  tae  problem  of  the  location  of  a 
moving  object  is  reduced  to  tae  integration  of  the  same  system  of 
differential  equations  (6.5),  waace  «x(t)  and  u^t)  are  considered  to 
be  known  functions  of  tiw  on  tae  oasis  of  formulae  (10.1),  (10.2) 
and  (10.6). 
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